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Polynomial  Distributed  Lags 


In  general,  a  distributed  lag  model  can  be  written 

p-1 

(1)     y^  =   2    P^x  T      ; 

^        T=0       ^ 


y,  and  x,  are  time  series  and  p   are  the  coefficients  of  the  lag  function.  Often 
the  number  of  periods,  p,  covered  by  the  lag  fionction  is  so  large  that  the  indiv- 
idual coefficients  3  cannot  be  estimated  with  sufficient  accuracy.  In  this  case 
we  usually  seek  to  estimate  the  coefficients  subject  to  some  restrictive  hypothesis 
a  priori .  The  best-known  hypothesis  is  that  of  Koyck  and  many  others,  requiring 
that  p  have  the  special  form 

(2)        p,    =    :^" 


A  great  deal  of  effort  has  been  devoted  to  the  econometric  aspects  of  estimating 
Koyck  distributed  lags,  without  achieving  any  general  agreement  about  how  to 
estimate  them.  The  purpose  of  this  paper  is  to  describe  an  alternative  lag  specific- 
ation which  is  both  more  flexible  and  easier  to  estimate.  Credit  for  the  discovery 
and  introduction  of  this  method  goes  to  Shirley  Almon  (1);  this  paper  merely 
restates  her  method  (with  some  modifications  due  to  Charles  Bischof f  (2 ) )  in  a 
somewhat  simpler  form,  and  goes  on  to  describe  the  computer  implementation  of  the 
method . 
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The  basic  hypothesis  of  the  method  is  that  the  lag  distribution  p  is  a 

T 

smooth  function  of  the  lag  T  .  Our  interpretation  of  the  notion  of  smoothness  is 
that  the  function  can  be  approximated  closely  by  a  polynomial  of  fairly  low  order; 
that  is,  we  suppose  that 

.2   .         .     r.     _N-1 


(3)      P^   =  a  +  a  T  Ha„T^  +  ...  +  ay 
T       12    3  N 


where  N  is  usually  less  than  6.  Polynomial  approximations  of  this  kind  have  a  long 
history  in  numerical  analysis,  but  the  methods  of  polynomial  interpolation  developed 
in  connection  with  them  do  not  appear,  in  retrospect,  to  have  any  usefulness  in 
approximating  distributed  lags.  It  is  Mrs.  Almon's  use  of  Lagrangian  interpolation 
polynomials  rather  than  ordinary  polynomials  that  makes  her  presentation  of  the 
method  somewhat  obscure.  The  simpler  approach  of  the  present  paper  is  formally 
equivalent  to  Mrs.  Almon's  approach. 

The  polynomial  approximation  gives  rise  to  a  straightforward  linear  esti- 
mation problem.  We  begin  by  substituting  equation  (3)  into  (l): 

(4)     y,   ^   ""i:   («i-^  V^  ••••*-V''"'K-t 

T=0 

p-1  p-1 

=  «1  (t=0  Vt  )   -^   ^  ^Iq     "Vt  )  ^  ••• 


By  defining  new  va.riables  z,  .  which  are  moving  averages  of  the  original  variables, 

t,  J, 


as  follows: 


(5)      z,  .   =    E   T   X 
"t,j 


p-1   j-1 

E    T    > 
T=0        t-T 


•we  have  a  linear  model  of  ordinary  form: 

All  estimation  methods  which  are  appropriate  for  linear  equations  are  available  to 
the  investigator  of  distributed  lags  if  the  method  of  polynomial  approximation  is 
used. 

In  practice,  the  method  proceeds  as  follows.  First,  a  polynomial  weighting 

matrix,  A,  is  generated;  it  is  convenient  to  normalize  it  so  that  the  lag  interval 

lies  between  — -r  and  —77  (instead  of  between  0  and  p-l)  —  this  is  done  by 
p+1      p+1  ^  ^      ^  ^ 

T  +  1 

taking  powers  of  -"rj-  instead  of  powers  of  T  .  The  matrix  A  has  N  columns  and 
p  rows  ;  each  column  gives  the  weights  applied  to  the  lagged  x's  in  generating  one 
z-variable  while  each  row  gives  the  weights  applied  to  the  estimates  of  a  in  cal- 
culating the  lag  function  P  . 

Next  the  z-variables  are  generated  from  A  and  x  using  the  product  relation 
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(7) 


[^,1  •••  ^,n] 


't^n-l'*"*'    t-pfl 


1       -^       (     —     )^ 
p+l        ^      p+1      ' 


1       ^  {^f 

p+1  ^   p+1    ' 


_2_ 
p+1 


^   p+1    ' 


(   -^   )  N-1 
^   p+1    '' 


(   ^   )  N-1 

^   p+1    ' 


N-1 


or. 


(8) 


x^A 


This  process  is  called  "scrambling"  at  MIT. 

The  next  step  is  to  obtain  estimates  a  of  the  coefficients  of  the  z-variables 
by  whatever  method  is  appropriate  for  the  stochastic  specification  chosen  for  the 
model.  Finally,  estimates  ^  of  the  lag  coefficients  can  be  obtained  from  a   using 
the  relation 


(9)      %       =       Act 


An  estimate  of  the  variance-covariance  matri 


X  V(^) 


can  be  calculated  as 


(10)  V(p)        =     A  V(a)  A^ 
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Two  additional  statistics  may  be  of  interest.  These  are  s,   the  sum  of  the 
lag  coefficients,  andya  ,   the  mean  lag.  If  u  denotes  a  vector  of  p  I's,  we  can 
calculate  s  from 


(11) 


A 

u'P 

A 

u'Aa 


and  its  variance  from 


(12)      V(s)   =   u'AV(^)A'u 


Second,  if  we  define  the  vector  v  by 


(13) 


V  = 


0 

1 


p-1 


then 


(U)     /I 


-  v'P 
s 


Its  asymptotic  variance  is 


(15)      V(ylL) 


(  Jv  -  ^2  u')AV(&)  A'  (i  V  -  J  V  -  \  u  ) 


TlB  process  of  calculating  the  lag  coefficients  and  these  statistics  is  called 
"unscrambling"  at  MIT. 
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We  turn  now  to  variations  of  this  basic  method;  these  take  the  form  of 
what  Bischoff  (2)  calls  zero  restrictions.  A  zero  restriction  is  used  to  impose 
a  priori  the  hypothesis  that  the  lag  distribution  approaches  zero  at  one  or  both 
ends.  If  the  method  of  polynomial  approximation  is  used  in  estimating  the  lag 
distribution,  zero  restrictions  are  imposed  by  limiting  the  components  z,  .  to 
those  corresponding  to  polynomials  which  meet  the  restrictions. 

Now  in  normalized  form,  the  basic  polynomial  lag  function  (3)  is 


(16)        P.    ^  «i    -   «2  i-~t-)    -   a^  ("irr-)'  ^ 

,J[!i_  ,  N-i 


If  a  zero  restriction  is  imposed  at  the  near  end,  formula  (l6)  is  modified  by 
eliminating  the  constant  term: 


,  T+1  sN-1 


This  form  of  the  distributed  lag  function  always  has  small  coefficients  for  the 
shortest  lags.  The  name  "zero  restriction"  is  derived  from  the  fact  that  if  a 
hypothetical  p   were  calculated  from  formula  (17),  it  would  be  zero  no  matter 
what  values  the  a-coefficients  had. 
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A  zero  restriction  is  imposed  at  the  far  end  in  a  similar  way.  Instead  of 
formula  (16),  we  use 


(18)    P 


a. 


+  ...  +    a 


T+l 

p+1 


+  a. 


-,  ^   T+l 
(111)2  _  

Vi       p+1 


N-1 


,I+lv  N-1 

Vl         P+1 


In  this  case,  a  hypothetical  p  is  always  zero,  so  that  the  lag  function  is  cons- 
trained  to  be  close  to  zero  for  the  longest  lags. 

Finally,  zero  restrictions  may  be  imposed  at  both  ends  by  dropping  the  first 
term  from  equation  (18): 


(19)         p 


a. 


(1+1)2 


T+l 
P+1 


+  a 


N-2 


xT+l,,N-l 

^p+1^ 


+     a. 


T+l 

P+1 


(1+1)3 

Vi 


T+l 
p+1 


Equation  (19)  is  the  form  which  Mrs.  Almon  proposed. 

Note  that  all  three  kinds  of  zero  restrictions  are  simply  modifications  of 
the  weighting  matrix  A  —  all  of  the  formulas  given  earlier  for  the  unscrambling 
phase  still  hold  for  the  modified  versions  of  A. 
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Programs  for  Polynomial  Distributed  Lags 

1.  AMAT 

This  routine  generates  the  A-matrix. 

Calling  sequence:      CALL  AMAT(NPER,   NDEG,   JZERO,A) 

NPER        Number  of  periodsj  p. 

NDEG         Number  of  terms  in  approximating  polynomial.  On  input  it 

should  be  the  highest  power  of  -^  plus  one.  On  return  it  will 
be  reduced  to  take  account  of  zero  restrictions  so  that  it  is 
equal  to  the  number  of  columns  in  A. 

JZERO     Zero  restriction  code; 

1  for  both. 

2  for  far  only. 

3  for  near  only. 
4-  for  neither. 

A        A-matrix,  packed  by  columns;  length  =  NPER^NDEG. 


2,  SCRAMB 

This  routine  generates  the  z-variables,  given  the  matrix  A  and  the  variable  x. 
Calling  sequence:  GALL  SCRAMB(NOBX,NPER,NDEG,A,X,Z) 

NOBX     Number  of  observations  in  X.  Number  of  observations  in  Z  is 
NOBX  -  NPER. 
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NPER  Number  of  periods^  p. 

NDEG  Degree  as  returned  by  AMAT. 

A  A-matrix. 

X  Input  vector;  length  =  NOBX. 

Z  Output  matrix.   Length  =   (NOBX  -  NPER>NDEG. 


3.  UNSCRM 


This  routine  calculates  various  statistics  which  are  useful  in  interpreting 
distributed  lag  estimates. 

Calling  sequence:  GALL  UNSCRM(NPER,NDEG,N0V,L0C,DMEAN,SMEAN,SUM,SSUM,B,D, 
V,VS,A,S,W) 
NPER     Number  of  periods,  p. 
NDEG     Degree  as  returned  by  AMAT. 
NOV      Total  nvimber  of  variables  in  regression  including  those  not 

associated  with  the  distributed  lag. 
LOG      Position  in  B  of  the  set  of  estimates  for  the  distributed  lag 

being  unscrambled. 
mmm  Mean  lag, 

S/EAK     Standard  error  of  DMEAN. 
SUM      Sum  of  lag  coefficients. 

SSM     Standard  error  of  sum  of  lag  coefficients. 
B        Vector  of  regression  estimates.  Length  =  NOV. 
T)  Vector  of  estimates  of  the  distributed  lag.  Length  =  NPER. 
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V        Estimate  of  variance-covariance  matrix  of  B.  Length  =  N0'\fe<-2  . 

VS       Variance-covariance  matrix  of  estimates  of  a  for  this  distributed 

lag,  extracted  from  V  by  the  program.  Length  =  NDEO«t2. 
A        Polynomial  weighting  matrix  A.  Length  =  NPER^NDEG. 
S        Vector  of  standard  errors  of  D.  Length  =  NPER. 
W        Scratch  vector.  Length  =  NPER. 
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